Suppression of molecular decay in ultracold gases without Fermi statistics 
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We study inelastic processes for ultracold three-body systems in which only one interaction is 
resonant. We have found that the decay rates for weakly bound molecules due to collisions with 
other atoms can be suppressed not only without fermionic statistics but also when bosonic statistics 
applies. In addition, we show that at ultracold temperatures three-body recombination involving a 
single resonant pair of atoms leads mainly to formation of weakly bound molecules which, in turn, 
are stable against decay. These results indicate that recombination in three-component atomic gases 
can be used as an efficient mechanism for molecular formation, allowing the achievement of high 
molecular densities. 

PACS numbers: 34.50.-s,34.50.Cx,67.85.-d,31.15.x 



In recent years, the efficiency of diatomic molecule for- 
mation in ultracold quantum gases and their stability 
once formed have become key ingredients for many exper- 
iments. For example, the remarkable stability of weakly 
bound molecules in ultracold gases of fermions in two 
different spin states [1] has greatly helped experimental 
studies in the BEC-BCS crossover regime 2] when tuning 
the two-body s-wave scattering length a from a > to 
a < near a Feshbach resonance. The long lifetimes re- 
cently observed for ultracold hcteronuclear molecules [3] 
might also pave the way for future studies on ultracold 
polar molecules. The stability of weakly bound molecules 
where at least one atom is a fermion has been inter- 
preted as a result of suppression due to the fermionic 
statistics governing both atom-molecule and molecule- 
molecule collisions d, . When bosonic statistics plays 
a role, however, weakly bound molecules tend to decay 
rapidly, limiting the experimental possibilities for studies 
near a Feshbach resonance [5, 6, 7]. Obviously, besides 
stability, efficient production of molecules is also of cru- 
cial experimental importance. 

In this Letter, we present a series of new predictions 
that opens up new possibilities for obtaining long molec- 
ular lifetimes and that also provide an alternative path 
for efficient molecular formation. Surprisingly, we have 
found that weakly bound molecules of two resonantly in- 
teracting atoms can still be stable against collisions with 
other atoms even when bosonic statistics applies. Specif- 
ically, inelastic collisions between these weakly bound 
molecules and any atom not resonant with either of the 
molecule's atoms is suppressed as a~^, therefore allowing 
long molecular lifetimes near a Feshbach resonance. This 
suppression in the absence of fermionic statistics has re- 
cently been observed in Ref. in an ultracold mixture 
of weakly bound *^Rb^°K molecules and ""^K atoms in 
an nonresonant state. Here, we show that the mecha- 
nism that determines the suppression can be traced 
to the Efimov physics governing the atom-molecule in- 
teractions [5|. In contrast to the cases with fermionic 
suppressioii(3, 0, our predicted suppression holds 



for any mass ratio between the collision partners. The 
sole requirement is that only the atoms bound in the 
molecule interact resonantly. 

We will also show that three-body recombination in 
such systems displays interesting and potentially useful 
behavior. In the one- and two-component bosonic gases 
and in bose-fermi mixtures studied to date, recombina- 
tion can lead to a substantial fraction of deeply bound 
molecules 0, S, @, giving them in the process enough 
kinetic energy to escape from typical traps. Moreover, 
any weakly bound molecule formed can eventually decay 
rapidly due to collisions with other atoms and molecules 
[10|. The situation is better for fermionic gas mixtures 
of two spin species: molecules are stable [j, Q and the 
recombination rate into weakly bound molecules (cx Ta^) 
dominates recombination into deeply bound molecules 
These two ingredients ensured the 
achievement of highly efficient molecular formation and 
long molecular lifetimes However, because recom- 

bination for fermionic systems scales as Ta^, efficient 
molecule formation is restricted to relatively high tem- 
peratures or large scattering lengths [Tlj . In this regime, 
though, collisions can potentially lead to high molecular 
dissociation rates (oc T^a^ 5]), limiting the molecular 
lifetimes. 

As we will show, a three-component gas of atoms 
with only a single resonant interaction retains the ad- 
vantages of these fermionic systems — efficient forma- 
tion of weakly bound molecules and long lifetimes — 
with the additional advantage of having a non-zero re- 
combination rate at T ^ 0, allowing the formation of 
molecules with much lower temperature which is impor- 
tant to explore the quantum degenerate regime. We 
have found that for these systems three-body recombi- 
nation into weakly bound molecules scales as a* at ultra- 
cold temperatures while recombination into deeply bound 
molecules scales as a^. This weaker a-dependence en- 
sures that near a Feshbach resonance recombination will 
lead mainly to formation of weakly bound molecules. As 
shown in Refs. III, these weakly bound molecules 
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can stay trapped due to the small energy released, min- 
imizing atomic and molecular losses. The most impor- 
tant aspect of molecular formation using this mechanism, 
however, is that even in the absence of fermionic statis- 
tics, the molecules formed are stable against decay since 
the decay rate scales as a^^ . This combination of higher 
formation rate and lower decay rate maximizes molecule 
production. 

We study these ultracold three-body collisions by solv- 
ing the Schrodinger equation in the adiabatic hyper- 
spherical representation. The three-body inelastic col- 
lision rates are determined by solving the hyperradial 
Schrodinger equation given by (in atomic units) 
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Here, fi is the three-body reduced mass, E is the total 
energy, F,^ the hyperradial wave function, and v a col- 
lective index that represents all quantum numbers neces- 
sary to label each channel. Since the hyperradius R gives 
the overall size of the system, this equation describes the 
collective radial motion under the influence of the effec- 
tive potential Wi, , with inelastic transitions driven by the 
nonadiabatic couplings Vi,i/' f5|. 

In the regime where the two-body interactions are res- 
onant, i.e., when \a\ 3> tq (where is the characteristic 
range of two-body interactions), the effective potentials 
Wi, become universal, and several analytical properties of 
ultracold three-body systems can be derived. In fact, in 
our framework [sj both the a and the E dependence of the 
three-body collision rates can be determined by simply 
identifying the attractive or repulsive character of W^{R) 
as determined by the influence of Efimov physics [al. Pre- 
vious work in ultracold three-body collisions 0)M S @| 
has analyzed the influence of Efimov physics for systems 
where at least two of the pairwise two-body interactions 
are resonant. In this paper, we will show that Efimov 
physics still has an important impact even when just one 
interaction is resonant. 

For systems with single resonant interaction, the effec- 
tive potentials are repulsive in the range rp ^ i? ^ |a| 
and can be conveniently parametrized by the coefficients 
Po and pi, 01 : 
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For a > 0, pq is associated with the atom-dimer channel 
and therefore plays an important role in atom-molecule 
processes. The effective potentials associated with p„ de- 
scribe three-body continuum channels. The coefficients 
Po and Pi, in Eq. ([2]) can be determined by assuming 
a zero-range model potential for the interatomic inter- 
actions (we have confirmed that finite-range two-body 
interactions give the same results) and writing the three- 
body wave function in terms of the Faddeev components 



[ill. Doing this, we have determined that po and Pi, are 

Po = J+l, p^ = J + 3,J + 5,... (3) 

where J is the total orbital angular momentum. Inter- 
estingly, in contrast to systems where two or three inter- 
actions are resonant, the strengths of the potentials pq 
and Pi, do not depend on either the mass ratio or the 
permutational symmetry of the resonant pair. 

Our derivation of these rates' scaling behavior rests on 
the observation that the rate limiting step is tunneling 
through the potential barriers represented by Eq. ^ . 
This, together with po and p^, from Eq. gives the scal- 
ing laws shown in Table [J for systems with one resonant 
pair of atoms. We show the results for XYZ systems 
(three distinguishable atoms) and for BBX systems (two 
identical bosons and a third distinguishable atom). For 
XY Z systems, we have assumed that the X-Y pair in- 
teracts resonantly and the Z atom is in an nonresonant 
state. For BBX, we assumed that only the intraspecies 
interaction is resonant and, therefore, the X atom is in an 
nonresonant state. Note that the scaling laws for both 
XY Z and BBX systems are exactly the same. In the 
table, we show the first three partial wave contributions 
to vibrational relaxation of weakly bound molecules, 
XY* + Z XY + Z (Frci); to three-body recombination 
into weakly bound molecules, X + Y + Z ^ XY* + Z 
(KJ)] to recombination into deeply bound molecules, 
X + Y + Z XY + Z (i^f ); and to dissociation of 
weakly bound molecules, XY* + Z ^ X + Y + Z {D3). 



TABLE I; Energy and scattering length dependence for three- 
body collision rates for systems writh a single resonant in- 
teraction. Boldface indicates the leading contribution at 
threshold (fc ^ where k = y/2^E for K3, and D3, and 
k — ^2pu(E — Ey) for Kci, —E^ is the molecular binding 
energy). The results are valid for all three-body systems 
with three distinguishable atoms (XYZ) or with two iden- 
tical bosons (BBX) with ubb resonant. 
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Figure [Ha) shows our numerical results for V^ci 
tained from a direct solution of Eq. ([T]) (see Ref. 
for details) for three distinguishable atoms, XYZ, for 
J'^ = 0+ (tt is the total parity). We have assumed a 
finite-range model two-body interaction in which the res- 
onant X-Y interaction is adjusted to support one weakly 
bound s-wave state whose position controls the scattering 
length, one deeply bound s-wave state, and one deeply 
bound p-wave state. Our model X-Z interaction sup- 
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ports only one s-wave state, and, for simplicity, we as- 
sumed the Y-Z interaction has no bound state. The fig- 
ure shows that when a/r^ > 10, our predicted scal- 
ing law is verified by the numerical results. Thus, we 
have clear evidence that vibrational relaxation is indeed 
suppressed for large a even in the absence of fermionic 
statistics. By extension, this result also support our ar- 
guments that molecular decay in BB + X collisions are 
also suppressed as even though bosonic statistics ap- 
ply. In a two-component gas of B and X atoms, however, 
long lifetimes for BB molecules can only be achieved if 
there are no free B atoms and with the additional con- 
dition of low molecular density, to prevent decay due to 
molecule-molecule collisions. 

Recombination for systems involving nonresonant 
atoms also presents some peculiar properties. As men- 
tioned previously, for a > 0, recombination into weakly 
bound molecules {K'^) scales as while recombination 
into deeply bound molecules {K^) scales as (see Ta- 
ble Therefore, recombination for large a leads mostly 
to the formation of weakly bound molecules. This dif- 
ferential scaling with a is in stark contrast to previously 
studied systems with two or three resonant interactions 
for which recombination is non-zero at T — 0, as is the 
case for one- and two-component bosonic systems and 
boson-fermion mixtures. In these systems, recombination 
into all two-body states scales as a^. In those systems, 
recombination into the weakly bound state is typically 
dominant as well, but only via the coefficient of . Sys- 
tems with only a single resonant interaction share this 
property in addition to the differential scaling with a. 
Figure [Ijb) shows our numerical calculations of recom- 
bination rates for the same system as in Fig. [TJa). The 
numerical results support all of the above predictions, in- 
cluding not only the smaller magnetude of relative to 
K'^ but also their scaling with a. Another consequence 
of this differential scaling is that for a < 0, i.e., in the 
absence of weakly bound molecules, recombination is pro- 
portional to |ap. So, recombination at a < for systems 
with a single resonant pair could be suppressed compared 
to a > for large a. 

The results obtained here for three-body processes in- 
volving a single resonant pair — summarized in Table |T] 
— combined with previous results [H] exhaust all pos- 
sible three-body processes relevant to three-component 
atomic gases. The competition between the various colli- 
sion processes in a multi-component gas is dictated fun- 
damentally by their energy and scattering length depen- 
dence, although it is also crucial to take into account the 
atomic and molecular densities. In fact, manipulating 
the atomic and molecular densities can make favorable a 
particularly desired coUisional behavior. 

For completeness, we show in Table HI] the three-body 
collision rates for all possible ultracold three-component 
atomic gases. For each mixture, we have listed the three- 
body processes according to which interatomic interac- 
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FIG. 1: (a) Vibrational relaxation rate from weakly bound 
XY* molecules and (b) three-body recombination rates into 
weakly, K'^ , and deeply, K^, bound molecules. 



tion is resonant by indicating the resonant pair of atoms, 
and we have only included the dominant contribution at 
threshold. Three-body processes involving a single reso- 
nant pair are indicated by "o" while all the others were 
obtained in previous work 0, 0, @]- Futher, we have 
2 < Po < 4 for Vtci and < po < 2 for depending 
on the mass ratio between the collision partners , and 
we have omitted modulation factors that produce the in- 
terference minima and resonant peaks related to Efimov 
physics [5]. 

We can see from the table that in all possible three- 
component mixtures weakly bound molecules can be ex- 
pected to be long-lived with relaxation scaling as a~^, 
provided that all free atoms are nonresonant. In fact, 
for gas mixtures where distinguishable fermions are res- 
onant, weakly bound molecules are expected to be long- 
lived irrespective of the identity of the free atoms. 

Given the suppression possible in atom-molecule col- 
lisions, it might be that molecule-molecule collisions 
are the dominant relaxation mechanism in a three- 
component gas. For BF molecules, that are themselves 
fermions, molecule-molecule collisions are suppressed at 
ultracold temperatures due to their p-wave character. 
For BiBj bosonic molecules, however, long lifetimes can 
only be reached for low molecular densities, in order 
to prevent molecule- molecule collisions. In contrast, 
for FiFj bosonic molecules, relaxation due to molecule- 
molecule collisions scales as a"'' with s > 1 Q for mass 
ratios mp./mp^ > 0.11603 (a condition satisfied for com- 
monly used alkali atoms) and the main molecular colli- 
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TABLE II: All three-body collision rates relevant for three-component atomic gases. For each mixture, the resonant pair of 
atoms are indicated by parentheses. The symbol "o" indicates collisions involving a single resonant interaction. 
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sional decay is due to atom-molecule collisions involving 
a single resonant interaction, scaling as a~^. 

From Tablellllwe see that the optimal three-component 
systems utilize molecular formation close to an F1-F2 Fes- 
hbach resonance. In this case, recombination involving 
two identical fermions is suppressed at ultracold tem- 
peratures while recombination in both B-F1-F2 and Fi- 
F2--F3 mixtures leads mainly to the formation of weakly 
bound F1F2 molecules. The most important property of 
this system is that after the weakly bound molecules are 
formed, molecular loss due to collisions with the remain- 
ing atoms is suppressed at least as a"^. Similar benefits 
can also be realized near a B-F resonance under certain 
circumstances. For instance, in B-F1-F2 mixtures, effi- 
cient B-Fi molecular formation could be achieved if the 
number of nonresonant F2 atoms is higher than the B and 
Fi numbers. If in the end all bosons are bound in BF 
molecules, then collisions with the remaining fermions 
is suppressed, ensuring the desired stability of weakly 
bound molecules. 

To summarize, we have shown that vibrational relax- 
ation of weakly bound molecules can be suppressed near 
a Feshbach resonance even in the absence of Fermi statis- 
tics. Remarkably, the root cause of this suppression can 
be traced to the same kind of universality that yields the 
Efimov effect. The systems that show this behavior can 
also exhibit differential scaling with a of the recombina- 
tion rate into weakly bound molecules that is finite at 
zero temperature. These two results combine to make 
ultracold three-component gases with a single resonant 
interaction an attractive alternative for efficiently pro- 
ducing stable, ultracold molecules. 
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